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Abstract 

The  follovlng  two-dimensional  problem  is   considered:      An   arbitrary- 
antenna  is   located  at  the  origin  which  laimches   surface  waves   on  a  reactive 
half-plane  y  =  0,    x  <  0.      It   is   assumed  that  the  other  half -plane,    y  =  0, 
X  >  0   is  occupied  by  a  perfect   conductor.      The   antenna  is   approximated  by 
a  distribution  of  multipoles  easily  found  from  the   antenna's   radiation 
pattern.      With  this   approximation,    simple   algebraic   expressions   are   ob- 
tained for  the   ratio  of  the  power  delivered  to  the   stirface  waves   as   com- 
pared to  the  power  in  the  radiated  cylindrical  field. 
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Figure  I 


2. 


2.   Formulation  and  Solution 

The  problem  (cf.  Fig.  1)  we  shall  solve  is  posed  by  the  three 
equations: 


(1)       (V^  +  k2)u(x,y)  =  -26(x)6(y), 


all  X,  y  g  0 


(2)       ^  +  Xu  =  0,  X  <  0,  0^  =  0 

(5)       1^=0,  x>0,  y  =  0 

vhere  u(x,y)  represents  the  transverse  magnetic  component  H  ;  H  =  H  =  0. 

For  convenience  k  has  a  small  positive  imaginary  part  k  =  |k|e  ,    and 

-icJt  * 

a  time  factor  e     is  surpressed.   Owing  to  the  normalization  chosen 

for  the  source  term,  the  free  space  Green's  function  of  (l)  is 


^        (i\  1    r   P^^^  +  i\/k  -V   |yl 

C        \/k' 


J?  ^ 


where  the   contour  C  is   taken  along  the   real- axis  in  the   v-plane  as   shown 
in  Figure   2.      Since   (k)    already  has  the   correct  singularity  at  the   origin, 
it   follows  that  if  we   seek  a  solution  of  the   form 


(5)  u(x,7)   =  ^   J    A+(v)     ^ 


.     /  2     2 
ivx  +  i  \/k  -V 


dv 


C  ^k^  -   V 


where 


(6)  A'*"(v)   =  0(1)    as    [v[  -  CO    ,    Im  V  ^   0 

then  equation   (l)   will   continue  to  be   satisfied.       At  the  boundary  we 

Recall  that  the  usual  source  strength  of  -4jt5(x)8(y)   leads  to  the 
Green's   function  jti  H   ^      (kp). 
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The  magnitude  of  S  is  greatly  exaggerated 
for  purposes  of  clarity 


Figure  2:    Complex  i/  plane 


assxane  that  we   can  differentiate   freely  under  the  integral  sign,    and 
then  if  A  (v)    is   an  analytic   function  in  the  upper  v-plane  the  boundary 
condition   (5)   will  be   satisfied  since  we   can  close   the   contour  C  for 
X  >  0.     Again,    it  is  sufficient  that  if  G"(v)   defined  by 


(7) 


G-(v)    =   A-"(v) 


1  + 


1  \/k     -    V 


is   an  analytic   function  in  the  lower   v-plane   and  O'(l)    as    |  v  |  -*  00    , 
Im  V   S  0  the  boundary  condition   (2)   will  be   satisfied.      Kay   ^ J ^    and 
Bazer  and  Karp    L  J    show  that  the   function  o(v)    can  be  expressed  as 


(8) 


a(v) 


_     ,^ 


a"(v) 


where   a   (v)    is   analytic   in  the  upper  v-plane   Im  v  g   im  k  and  a~(v)    is 
analytic  for  Im  v  S  Im  k.      It  follows  that  equation   (7)    can  be  written  as 


(9) 


G-(v)    a"(v)   =   a"(v)    a"(v). 


It  can  be  shown  that  \a   (v)|  and  \a   (v)|  are  bounded  above  and  below  by 
positive  constants  in  their  respective  domains  of  analyticity.   Owing 
to  the  postulated  behavior  of  g"(v)  and  A  {v) ,    the  left  side  of  (9)  is 
boxinded  and  analytic  for  Im  v  ^  0  and  the  right  side  is  bounded  and 
analytic  for  Im  v  a  0  implying  that  each  side  is  the  analytic  continua- 
tion of  the  other.   If  we  invoke  the  Liouville  theorem  we  discover  that  each 
side  of  (9)  is  equal  to  some  constant  K  proportional  to  the  strength  of 
the  source,  whence 


For  example  reference  [2] 


4. 


g"(v)    a"(v)   =  a'^(v)    a^(v)   =  +  K 


(10) 


A"(v) 


L  G"(v)   = 


a^(v) 


a"(v) 


Observe  that   (lO)   yields  kernels   A  (v)    and  G   (v)   having  the   required  asymp- 
totic behavior  for  large    |v|. 

The  integral  representation  of  the   solution   (5)    is   then 


(11)  u(x,y)  .  ^   /   - 


/~2     2 
1  e^^^  +  i  s/k   -V     y 


a'(v) 


v/£^^ 


dv 


or  alternately  by  {8) 


(12)  u(x,y; 


K        r    1 

2«i   J     .-rvv 


/2     2 
^ivx  +  i  \/k  -V     y 


dv 


o-i^)oiv)  j^T7 


From  (11)  we  observe  that  the  poles  of  the  integrand  are  the  zeros  of 
a  (v).   Since  ]a  (v)|  is  bo\mded  above  and  belov  in  the  upper  half  plane^ 
these  can  only  arise  for  Im  v  ^  Im  k.   For  reference  we  locate  the  poles 
of  the  Integrand  of  (ll) . 


We    calculate   that 


•     A  2  2 

1  Vk     -   V 


a  (v)        a"(v)a(v)  a"(v)        i  X^T^  +  x 

,  i  /k^  -  v^  (X  -  1  /k^  -  v^) 


^^)         (.A^7^.v)(yx2^k2.v) 


The  only  pole  that   can  arise  la  in  the  lower  v-plaae   at  v     given  by 


(13)       v^^  -J^TIF . 


The  corresponding  residue  of  (ll)  is  found  to  be 


(lU)  K >^      g-i  /k^  +  X^  X  -  Xy 

a"(-\/k^  +  X^)  /x^  +  k^ 


which  represents  the  suxface  wave  contribution. 

3.    The  Far  Field 

Introduce  the  transformations 

(1)  V  =  k  cos  Cp 

and 

(2)  X  =  p  cos  i3,  y  =  p  sin  * 
80  that  (2.11)  becomes 


6. 


f     i 

-L  ■*"(k  c( 


/'z^  ^n    o^  K        r  1  ikP   co3(cp-'a)   ^ 


^p  (k    COB    cp) 


in  which  the   contour  C  has  been  transformed  into  the  path  P   (Fig.    3)- 
If  ve  transform  the   alternate  representation   (2,12)    to  the  cp-plaae,   we 
obtain 


r.        n      .1,   „.-,  „  „ikP   cos(cp-«) 
(4)  u{p,i)   =   -^J        iksmcpe 


o   (k   cos  cp)(lk  sin  cp  +  ^) 


from  which  it  is   seen  that   the  pole   at    v     which  generates  the  surface 


(5)  sin  cpQ  =   -  ^  or  cos  cp^  =   -  /i  +    (\/k)^   . 

An  asymptotic  development  of  either  (5«5)  or  its  equivalent  (5«^)  is  readily 
obtained  by  deforming  the  path  P  into  the  steepest  descents  contour  S(  i3) . 
For  a  given  ^  this  contour  is  defined  as  the  locus  of  points  in  the  cp- 
plane  which  satisfy  the  equation 

2 

(6)  cos(cp  -  ■«)=>'  1  +  i  s 

with  s   as   a  real  parameter.      It  is   easily  established  that  S('6)    origi- 
nates  at  Re  cp  =   '6  -  «/2,    Im  cp  =  <»,    passes   through  the   saddle  point  cp  =   i8 
on  the   real  axis,    and  terminates   at  infinity  Re  cp  =  ■e  +  ■a/2,    Im  cp  =   -  ". 
Inspection  of  Figure   5  reveals  that  for  0  §  ■6  <  rt/2  we   can  deform  C 
freely  to  S('fl)   without  picking  up  the   residue   gt  cp  ,    and  for  some  t  in 


6a. 


Steepest    descents 
contour    S(^)   for 


Figure  3:     Connplex   (/>  plane 


7. 


the  Interval  rt/2  <  -3  S  n  we  pick  up  the   residue.      With  the   aid  of  (5) 
and   (6)    it   is   established  by  comparing  the  real  parts   of   (5.6)   that  this 
critical  angle  ^     is   given  by 


(7) 


x/l  +   (x/k)' 


where  jt/2  S  ■a     ^  it   so  that  we  only  obtain  the  surface  wave  on  the  left 
'  c 

for  X  <  0. 

For  kP  »  1,   we  then  obtain  the   far-field  behavior 


(8)  lim       uCp^'O)   =  K 

p    -♦    00 

where 

(9)  m=     J^+ 


-imx-Xy 


^ik(P  +  TtA) 


a~(-m)  m 


U(^-^^)   +  — 


a  (kcos^)  y2itk 


*  e,[(.p)-H 


^^ 


and  U(*)    is  the  Heaviside  unit  function.      The   result   (3-8)   should  be 
modified  somewhat,    if   X    -^0,    to   allow  for  a  pole  in  the  vicinity  of  the 
saddle  point  for  ^  =  «.      This   detail  is   fully  discussed  in    [2]    and    \^\ 
where  the  procedure   for  Isolating  the  pole   is  explained. 

Apart  from  terms  ©•[(kp)'^''^  j,    the   radiation  field  contains   a  cylin- 
drical wave  p(^)e^^'^/7kP  "^"tli  a  pattern  function  -gi^)   proportional     to 


(10) 


a~(-  k  cos  •6), 


a   (k  cos  iS) 
and  in  addition  a  surface  wave  field  whose  amplitude  is  proportional  to 

(11)      -^  . 
0   (-■) 


Since  a(  v)  Is  even  in  v,  we  have  0^{v)   =  l/a"(-v),  [2]. 


8. 


For  each  field  the  power  carried  off  to  infinity  varies  as  the  absolute 
square  of  these  coefficients.   As  Kay  and  Lurye  have  observed,  it  is  an 
extraordinary  occurrence  in  that  the  absolute  squares  of  (5-10)  and  (5-11) 
axe   simple  elementary  functions.   From  Kay's  paper,  equation  (66),  we  ob- 
tain in  our  notation: 

(12)  |0"(-k  cos  ^)|^  =     ''"^  ^  "^  ^        -  OS^^n 


cos   ■8 


+  yi  +  (x/k)2 


72  2 

(15)  |a-(-m)r  =    |a-(-yk''+   X^)  |      = 


2j^V? 


2m 


Note  that  our  X  corresponds  to  Kay's   a. 

k.        The  Energy  Distribution 

In  this  section  we  shall  obtain  an  explicit  expression  for  the  ratio 
of  the  power  delivered  to  the  surface  wave  field  as  compared  to  the  cylin- 
drical radiation  field.      From  (3.8)   we   see  that  the   surface  wave  is 


(1)  u_^(^.y) 


s\irf'    '■"         ^      -f       ^ 
m  a   (,-m; 

The  Poynting  flux  [T  transported  by  this  field  is  readily  calculated 
we  recall  that  u(x,y)  represents  the  transverse  magnetic  component  H  . 


9. 


n  surf  =    /     fy-^Ijsurf  ^^  =  it^    f       ^  ^^^^^  ^^ 


X         I       K 


,2 


a   (-m) 


If  we   use  Kay's   result    (3.15)    then 

(2)  n  ^^^^  =     AkL =   tsnhi^  sinh  A) 

UK  (Jl   +  A   +  ij 

vhere  A  =  x/k. 

In  the  same  fashion  we  can  compute  the  radiated  power  of  the  cylin- 
drical field  by  using  (2.12) 


a   (k   cos  •6  )  y2rtkP 


cyi     J       V    .  ^^^ 


-  f  h<^ 

O 


pd^ 


1^1      r  d^ ^  jfL    r        cos  ^  +  1 

2rtU)e     /       I    +,,  „^|2     ~   2rtU€    J  / ; 

o      1^   ^^  '^^^  ^^'  o     cos  ^  +yi+   (x/k)' 


The  integral  that  appears  in  this  expression  is  evaluated  in  Appen- 
dix A;  we  find 


w        n  =yi  =  Sir 


2  + 


-  Vl  +  A^  + 


A(A  -  yi  +  a" 


The   ratio  J]         f  /  TI  represents  the  efficiency  of  launching  the   si 

face  wave.      With  the   aid  of  {2.k-)    and   (^.^)   we  obtain 


(5) 


n 


surf 


ll  cyl    (1  +  Vl  + 


- 

a(a  -    Vl  +  a  ) 

2A(A   - 

•y/l  +  A     +A-     v/l+A     +1 

For  A  »  1,  the  bracketed  exrpession  behaves  like  2A  so  that 


(6) 


'--    Ilcyl 


Consequently,  even  without  using  a  directive  soiirce  the  major  portion  of 
the  radiated  energy  is  apportioned  to  the  surface  wave  as  we  increase  the 
surface  reactance  A. 


It  is  perhaps  of  greater  use  to  conrpaxe  the  surface  wave  energy  to 
the  total  radiated  energy 


11  total       11 


total       11  cyl       li  surf 


If  we  use   (2.^4-)    and   (^.^)    a  little  manipulation  reveals 


lOa. 


A 

li  surf  /  li  tot 

neyi/ntot 

li  cyl  /  li  surf 

0.00 

0.0000 

1.0000 

00 

0.02 

0.0198 

0.9802 

49,5050 

0.04 

0.0592 
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24.5100 

0.06 

0.0582 

0.9418 

16.1816 

0.08 

0.0768 

0.9252 

12.0200 

0.10 

0.0950 

0.9050 

9.5240 

0.12 

0.1128 

0.8872 

7.8652 

0.l4 

0.1502 

0.8698 

6.6777 

0.16 

0.1^75 

0.8527 

5.7897 

0.18 

0.1659 

0.8561 

5.1002 

0.20 

0.1802 

0.8198 

4.5495 

0.22 

0.1961 

0.8059 

4.0998 

0.24 

0.2116 

0.7884 

5.7258 

0.26 

0.2268 

0.7752 

5.4101 

0.28 

0.2415 

0.7585 

5.l4oi 

0.50 

0.2560 

0.7440 

2.9067 

0.52 

0.2700 

0.7500 

2.7051 

0.5^ 

0.2858 

0.7162 

2.5259 

0.56 

0.2972 

0.7028 

2.5650 

0.58 

0 . 5102 

0.6898 

2.2254 

o.J+o 
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0.^2 

0.5554 
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O.kk 
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0.6525 

1.8779 

o.h6 

0.5592 

0.6407 

1.7854 

o.k& 

0.5708 

0.6292 

1.6971 

0.50 

0.5820 

0.6180 

1.6180 

0.60 

0.4558 

0.5662 

1.5052 

0.70 

0.4795 

0.5207 

1.0862 

0.80 

0.5194 

0.4806 

0.9254 

0.90 

0.5546 

0.4454 

0.8050 

1.00 

0.5858 

0.4i42 

0.7071 

1.25 

0.6492 

0.5508 

0.5405 

1.50 

0.6972 

0.5028 

0.4545 

1.75 

0.7544 

0.2656 

0.5616 

2.00 

0.7659 

0.2561 

O.509L 

2.25 

0.7878 

0.2122 

0.2694 

2.50 

0.8074 

0.1926 

0.2585 

2.75 

0.8258 

0.1762 

0.2158 

5.00 

0.8577 

0.1625 

0.1957 

4.00 

0.8769 

0.1251 

0.l4o4 

5.00 

0 . 9010 

0.0990 

0.1099 

6.00 

0.9172 

0.0828 

0.0902 

8.00 

0.9577 

0.0625 

0.0664 

10.00 

0.9501 

0.0499 

0.0525 
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1 1^  total  u)€ 


A         .    A   -   M  +    1 
M  +   1        2A(A   -  M) 


M  =    v/l  +  a' 


(8) 

11  surf 
11  total 

A 

A             A   -   M  +   1 

M  +  1    '    2A(A   -  M)    '    ^ 

1  +     v/l  +  A^ 

which  is  plotted  in  Flgiire  h   as   a  function  A. 

5.        Arbitrary  Excitation 

In  actual  practice   one  vould  not  excite    a  surface  vave   structure  by 
a  line  source.      We  would  use   a  directive   antenna  A  such  els   an  exponential 
horn  to   focus  the  maximum   energy  into  the   surface  wave.      How  can  we   obtain 
the   excitation  produced  by  such  a  radiator?      If  we   replace   the   surface  wave 
structure  by  a  metallic  plane  then  the  antenna  A  will  produce  some   radiation 
pattern  j){-b)   which  will   characterize  the   antenna.      The  pattern  p(6)    is  known 
for  0  S  •a   S  jt    so  that  it   is    a  simple  matter    [Appendix  b]    to   approximate   it 
arbitra.rily  by  a  cosine  polynomial 


N 


(1)  p(^)    ~  v/^) 


a   (k.   cos   -Q) 


since  the   cosine  functions   cos  U'fl   or  cos  ■&   are   complete   on  the  interval 
0  S  ■e  S  rt.      In  particular  the   sine  functions   sin  U'd   are   expressible   as 
infinite  series  of  cosine  functions   for  0  g  ^  S  it.      We  observe  then  that 


11a. 


ay 


N      a 


(n) 


V^^  .  k'u  =-2S   -f   8     (x)    S(y) 

n  =  0      i 


iiiiiiiiiiiiniiiiiililllll 


ay 


+   Xu  =  0 


-2E^S^"^(x)S(y) 


n  =  0 


au 
ay 


=  0 


Figure  5:    Approximation  to  on  arbitrary 
surface  wave  launcher 


12. 


•     A  2  "      2 
+  i  Vk 


1         r      X~  n     e^^^  +  1  Vk      -    V     y 


is   a  ftmctlon  vhlch  possesses  the   approximate   far- field  "behavior  (l).      It 
is   also   clear  that   (2)    satisfies  the  wave  equation   (cf.   Figure  5) 


N 
"     %     (n).    .     .    . 

y  s  0 


IN 

(3)  (V^  +  k2)u(x,y)   =   -   2     21     ^  6^''^(x)5(y), 

n=o     i 

where  5^^'^(x)   =  S^/Sx^     5(x)    and  the  boundary  condition 
(k)  |H.o  y=0. 

Equation   (3)    reveals  that  what  we  have  done  is   to  replace  the   antenna  by 
an  approximate   equivalent-   a  distribution  of  nrultipoles.     We  shall  assume  that 
the  a     are   real  which  corresponds  to  feeding  the  multipoles  in  phase   inter- 
vals of  rt/2.      One  need  only  solve  equation   (3)    subject  to  the  two  part 
boundary  condition 

(a)        ^=  0  y=  0,    X  >0 


(5)  \ 


(b)        ^  +  Xu  =  0  y  =  0,    X  <  0 


to  obtain  the  effect  of  introducing  the  surface  wave  structure  in  the  left 
half  plane. 


15. 


Our  procedvire   is   to   choose 


+   i\A 


2  2 


(6)  u(x,y)=^     /      A-'(v)E(v)      ^ — IZIZ: 

J  r2        2 


vhere  E(v)  is  a  known  entire  function  of  v,    I.e.,    a  Chebyshev  polynomial 

apprDximation  defined  by 
N 
(   E(v)  =  21  ^n^"" 


(7)     < 


E(k  cos  ^)  =  p^(i3) 


If  A  (v)   is   an  analytic  function  of  v     for  Im  v  S   0,    the  boundajry  condi- 
tion  (5a)   will  be   satisfied.      In  addition,    if 

(8)  A-'(v)=Otl).  |v|-», 

then   (6.5)   will  have  the   singularity  at  the  origin  demanded  by  (5.5). 
For  X  <  0,    y  =  0  we   observe   that  if 


(9)  G-(v)   =  A-'(v)E(v; 


is  an  analytic  function  of  algebraic  order  in  the  lower  v-plane,  (5b)  will 
be  satisfied.   If  we  assume  (8)  then  we  see  by  (9)  that 


(10)      G"(v)  =  0(v") 


li^. 


Use  of  (2.8)  yields  by  (9) 

(11)      G-(v)  a-(v)  =  A+(v)  E(v)  a^{v) . 

Once  again  we   repeat  the   airgument  that  each  side  of  (ll)    is   analytic  in 
a  half  plane  with  a  common  region  (at  least  the   real  axis)    of  analyticity 
and  hence  is  the  analy-tic   continuation  of  the  other.      Owing  to  the   assiuned 
growths ;,   we  obtain  by  an  appeal  to  the  Liouville  theorem  that  each  side   is 
equal  to  a  polynomial  P   (v)    of  order  N,    hence 


a'^(v)   E(v)    a"^(v)   =  P^v) 


+  ...      +...        fh) 


(12)  A"(v)    a'(v) 


Since  E(v)  is  a  polynomial  with  real  coefficients  its  zeros  are  either  on 

the  real  axis  or  conjugate  pairs  on  the  imaginzary  axis.   Unless 

P  (v)  =  K  E(v)  where  K  is  some  constant,  the  right  side  of  (12)  will  contra- 


(15)     A"(v)  = 


a'"(v) 


m  u(x,y)  =  ^  /  ^ 


^iv. 

:+iyk^-v^; 

Y 

y.^  -  v^ 

and  again  the  constant  K  is  proportional  to  the  strength  of  the  source. 
A  straightforward  application  of  the  method  of  saddle  points  yields 


15. 


(15)  llm       u(p,,«)   =  u^^^^  u(^   -  ^^)  +  u         +  Cr[(kp)"5/^J 

P    -►    00  "^ 


,^r\  ,        \       XK  E(-m)        -i  yk^+  X^  x-Xy 

a   (-m) 


and  the   cylindrical  field  u  is 

KP,(.)  e^'^'^^^ 

^y^       o  (k  COS  «)         y^tkp 


The  power  pattern  P(^)  of  the  cylindrical  field  is  distributed  as  a 
function  of  ^  by 


(18)     P(.)=  |p,(.)l^  -    ^°^^^l 


+  v/l  +  (X/k)^ 


It  is  interesting  to  observe  that  this  expression  is  |pj^(i8)|   modxxlated 
by  the  prevlcms  pover  pattern  function  of  the  line  so\xree. 

The  Poynting  fliix  j]  transported  by  each  field  is  readily  calculated 

I  1  surf   ico€  /  \_0y.        J  ^^^^  ''    we   m  +  k  '  '    ' 


From  (5-7)  ve  have  E(-m)  =  p^(cos~  ^)  =  p^(cos'  ^1  +  A  )  so  that  the 
amplitude  of  the  surface  wave  is  proportional  to  the  reactive  power 


16. 


-1    J 2 

delivered  by  the   antenna  to  the   complex  angle  •fl  =  cos      ^Jl  +  A   , 


The   cylindrical  radiated  pover  becomes 

1    f  1^1'  r 

IT  =    _i-        /     u    U*    d^    =    ^r / 

II  cyl        ia)€    J       P  23tto€    J 


vM 


(k  cos  ^; 


d^ 


|Kl 


2        :t 


(20) 


Ilcyl  =  S^/     IV^)|^ 


cos   ^ 


+  \/l  +  a' 


d-d 


and  then 


In  Appendix  A  ve   evaluate   Integrals   of  the   form 
I     =      /      cos   i 


cos  ^  +  1 


cos  i5 


+  v/l  +  a' 


vhlch  vould  materially  ease  the  labor  in  computing  (21)  for  any  specified 
radiation  pattern. 
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Appendlx  A 


Evaluation  of  the  integral  I  where  M  =    /l  +   (-) 

r      cos  «  +  1    ,,        1      r      cos  ^  +  1    ,,  l«l    >.  1 

^o  =    J        cos  ^  +  M  ^^=2    J        cos  ^  +  M  ^*   '  1^1    ^  ^- 

o  o 

This   Integral   can  be    evaluated  very  simply  "by  residues;    define 

i^ 
z  =   e 

vhence 

dz  =  i   z  d« 

cos  «  =  -(z  +  -) 

so  that  I     becomes 


iz  =1 


The  denominator  can  be  readily  factored 


(z^  +  2Mz  +  1)  =  (z  -  a)(z  -  p) 


vhere 

a  =  -  M  +  VM^  -  1  ,      3  =  -  M  -  yM"" 


and  since    \k\    >  1,    clearly    |a|    <     1,    \^\    >  1.      If  we   add  the  residues   at 
z  =  0  and  z  =  a,   ve  obtain 


18. 


I    g^  +   2a-  +  1  ^    l"1 


Ln  terms   of  A  =  X/k 


A  -     \/l  +  A     +1 


A(A   -     Jl  +  A^) 


It  is   a  simple  matter  to  note  that  I     -  0   as  A  -  <»,    and  I^  "*  n   as  A  -*  0. 
If  the  Integral  is   of  the   form 


r  °-Q      cos  ■6   +  1 

\  =    J      '^^^  ^      cos  ^  +  M 


d'O 


the  same  procedure  yields 


^  "  2«i  gn   y    z^  +  2MZ  +  1 


|z|=l 


Observe  that 


i(z  +  i)      =   1  +  i^ 
z^ 


2  2        2 

z^  z' 


^d  in  general  for  n  odd  the  coefficient  of  l/z  vanishes  for  this  factor; 


19. 


for  n  even  the  coefficient  of  l/z  is 


c  1^        r(n  +  1) 


Hence  for  n  odd,  the  only  pole  of  the  integrand  is  at  z  =  a,  we  find 


yi  + a' 


rA^  -  2A  Vl  +  A^  +  l\       .  . 


2,n  +  2s 


-yi.A^  y    '-(^^i^) 


n  odd 


20. 


Appendix  B 

Firsts  observe  that  the  cos  n^  functions  form  a  complete  set  on  the 
interval  0  ^  -d   S  jt.  Expand  p(i8)  as  a  Fourier  cosine  series 


(1)        p(«: 


Sa, 


in  the  interval  0  ^  i8  fe  jt  either  from  an  analytic  knowledge  of  p('6)  or 
nvunerically  l  -'  from  a  measured  radiation  pattern.   Tr\mcate  the  series 
at  some  n  =  N  and  express  each  cos  mS  as  a  polynomial  of  order  n  in 

cos  i,    e.g. 

4  2 

cos  J+'fl  =  8  cos  ■6-8  cos  •S  +  1. 

Such  a  collection  of  terms   in   (l)   will  yield 


(2)  p(«)   =     2     %  ^°^''^ 


or  as  a  function  of  v  =  k  cos  6 


(5)       Ejv) 


2L.  -n  ^k^ 


n=o 


The  expression  (3)  corresponds  to  a  Chebyshev  approximation  L J  of 
E  (v)  in  the  interval  |^|  S  1  by  a  polynomial  of  order  N. 
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The  Olenn  L.  Martin  Company 
Baltimore  3,    Maryland 
Attni      Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
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3000  Arlington  Boulevard 
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Wayland  Laboratory,    State  Road 
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Wayland  Laboratory 
Wayland,    Mass. 
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Farmingdale,   Long  Island,   N.  Y. 
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Research  Laboratory  of  Electronics 
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Cambridge  39,  Massachusetts 
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University  of  Michigan 
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Attn:  J.  A.  Boyd,  Supervisor 
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Engineering  Research  Institute 
Radiation  Laboratory 
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Engineering  Research  Institute 
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P621  Oeorgia  Avenue 
Silver  Spring,  Maryland 
ATTN:  Dr.  B.  S.  oiurary 
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Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 


Nuclear  Developnont  Associates,  Inc 

5  New  Street 

White  Plains,  New  York 

ATTN:  Library 

California  Institute  of  Technology 
Electrical  Engineering 
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P.  0.  Box  73 

Lexington  73,  Massachusetts 

ATTN:  Dr.  Shou  Chin  Wang,  Room  C-351 
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Laurence  0.  Hanscom  Field 

Bedford,  Mass. 
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University  of  Maryland 
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Washington  University 
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San  Diego  3,  California 
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Electrical  Engineering 
California  Institute  of  Technology 
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•TTM:  Dr Georges  G.  Weill 

Naval  Research  Laboratory 

Washington  25,  D.  C. 

ATTN:  Henry  J.  Passerinl,  Code  527eA 

Dr.  George  Kear 
5  Culver  Court 
Orlnda,  California 

Brooklyn  Polytechnic 

8?  Livingston  Street 

Brooklyn,  New  York 

ATTN:   Dr.  Nathan  Marcuvlti 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
fl?  Livingston  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Jerry  Shraoys 

Department  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexico 
ATTN:  Dr.  I.  "Colodner 

Mathematics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTN:  Dr.  Harry  Hoehstadt 

Ballistics  Research  Laboratory 
Aberdeen  Provl'-^g  Grounds 
Aberdeen,  Maryland 
ATTN:  Dr.  Pullen  Keats 

Dr.  Lester  Kraus 
Ii935  Whitehaven  Way 
San  Diego,  California 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,  Minnesota 
Attn:  Dean  Athelston  Spllhaus 


Naval  Research  Laboratories 

Washington  25,  D.  C. 

Attn:  W.  S.  Ament,  Code  5271 

Naval  Research  Laboratory 
Washington  25.  D.  C. 
Attn:  Dr.  Leslie  G.  McCracken,  Jr. 
Code  3933A 

Office  of  Naval  Research 
Department  of  the  Navy 
Attn:  Geophysics  Branch,  Code  Iil6 
Washington  25,  D.  C. 
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Signal  Plans  and  Operations  Division 
Attn:  SIOOL-2.  Room  20 
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Defence  Research  Member 
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Attn:  Technical  Reports  Library 
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U.  S.  Weather  Bureau 
U.  S.  Department  of  Commerce 
Washlnnton  25,  D.  C. 
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